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We present a first-principles study of bare and hydrogen passivated armchair nanoribbons of the 
puckered single layer honeycomb structures of silicon and germanium. Our study includes optimiza- 
tion of atomic structure, stability analysis based on the calculation of phonon dispersions, electronic 
structure and the variation of band gap with the width of the ribbon. The band gaps of silicon 
and germanium nanoribbons exhibit family behavior similar to those of graphene nanoribbons. The 
edges of bare nanoribbons are sharply reconstructed, which can be eliminated by the hydrogen 
termination of dangling bonds at the edges. Periodic modulation of the nanoribbon width results 
in a superlattice structure which can act as a multiple quantum wells. Specific electronic states 
are confined in these wells. Confinement trends are qualitatively explained by including the effects 
of the interface. In order to investigate wide and long superlattice structures we also performed 
empirical tight binding calculations with parameters determined from ah initio calculations. 

PACS numbers: 73.22.-f, 63.22.-m, 68.65.-k 



I. INTRODUCTION 

Graphene, a two-dimensional (2D) honeycomb struc- 
ture of single layer of carbon atoms, have been synthe- 
sized and established as a material with wide range of 
unusual properties It is a semimetal with a band 
profile having linear dispersion near the Fermi level, 4 
which attributes to its electrons a massless Dirac Fermion 
behavior. Quasi one-dimensional (ID) derivatives of 
graphene, called graphene nanoribbons, were also pro- 
duced recently!*™ Graphene nanoribbons are semicon- 
ductors having interesting electronic properties depend- 
ing on their geometry. These properties can be used 
to fabricate nanodevices as field effect transistors, spin 
valves, multiple quantum wells and etc.^UJ 

2D honeycomb struc tures and ID nanotubes of Si and 
Ge were studied ear lierP^E^ Stringent stability tests have 
recently shown that, 2D honeycomb structures of Si and 
Ge can be found stable in a slightly buckled geometry.^ 
These structures have similar properties as graphene and 
thus carry the potential of being used in the similar 
applications.^ Compared to graphene, the interatomic 
distance is larger in Si and Ge, so the diminished tt — ir 
overlaps cannot maintain the planar stability anymore. 
Eventually, the sp 2 hybrid orbitals are slightly dehy- 
bridized to form sp 3 -like orbitals, which in turn results 
in a puckered structured Freestanding graphene sheets 
and nanoribbons can be produced spontaneously, but it is 
not the case for 2D honeycomb structures of Si (silicene) 
and Ge. However, there are plenty experimental work on 
growth of Si nanoribbons especially on Ag surface 
These highly metallic nanoribbons are formed by self- 
organization and have straight, atomically perfect and 
massively parallel structures!^! The electronic structure 
of Si nanoribbons on Ag surface was also investigated 
theoretically!^ 

This paper reveals the atomic and electronic structure 
of armchair Si and Ge nanoribbons. Stability analy- 



sis, based on the calculation of the phonon dispersions 
via the force constant method, was performed. Energy 
band structure calculated by first-principles density func- 
tional theory (DFT) were used to generate the parame- 
ters of the tight-binding model. It is found that armchair 
nanoribbons of Si and Ge are stable and their band gap 
vary with their width displaying a family behavior! 20 * 21 * 
Formation of multiple quantum well structure in super- 
lattices consisting of periodically repeated junctions of 
nanoribbons having different widths was also investigated 
in detail. A ID model was proposed to understand the 
effect of the interface in superlattice structures. It is 
found that specific electronic states are confined in these 
superlattice structures. The interface effects, which can 
explain unexpected confinements, are revealed. 



II. METHODS 

We have performed first-principles plane wave calcula- 
tions within Local Density Approximation (LDA)P us- 
ing projector augmented wave (PAW) potential^. All 
structures are treated within supercell geometry using 
the periodic boundary conditions. A plane- wave basis set 
with kinetic energy cutoff of 300 eV is used. In the self- 
consistent potential and total energy calculations, the 
Brillouin zone (BZ) is sampled by (15 x 1 x 1) special 
k-points. All atomic positions and lattice constants are 
optimized by minimization of the total energy and atomic 
forces. The vacuum separation between the nanoribbons 
in the adjacent unit cells is taken to be at least 10 A. 
The convergence for energy is chosen as 10 -5 eV be- 
tween two steps, and the maximum Hellmann-Feynman 
forces acting on each atom is less than 0.02 eV/A upon 
ionic relaxation. Numerical plane wave calculations have 
been performed by using VASP package! 24 * 25 * Phonon dis- 
persions were obtained using the force constant method 
with forces calculated in a (5x1x1) supercellPSEBCou- 
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FIG. 1. (a) Atomic structure of fully relaxed bare armchair Si 
nanoribbon of width n — 9 (ASiNR-9). (b) Atomic structure 
and bond length distribution of hydrogen saturated ASiNR-9. 
Similar pattern is observed in hydrogen saturated AGeNR-9. 
The primitive unit cell of ASiNR-9 is shaded. The zigzag 
chain of Si atoms perpendicular to the nanoribbon axis is 
delineated by the dashed lines. 



pling parameters of empirical tight binding calculations 
are determined from the ab initio results and are used to 
treat wide and long super lattice structures, comprising 
as many as 3600 Si atoms. 
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III. ATOMIC STRUCTURE AND STABILITY 

We first investigate the atomic structure of bare and 
hydrogen saturated armchair nanoribbons of Si and Ge. 
The ideal honeycomb structure is cut parallel to the near- 
est neighbor bonds to form an ideal bare nanoribbon with 
a certain width. Armchair nanoribbons are classified by 
counting the number, n, of Si (or Ge) atoms forming a 
zigzag chain perpendicular to the cut direction. Accord- 
ingly, there are 2n Si (or Ge) atoms in the primitive unit 
cell of an ideal bare armchair nanoribbon. This structure 
is treated by a supercell having periodic boundary condi- 
tion in cut direction and a vacuum spacing in other direc- 
tions. To lift the constraints imposed by (lxl) unit cell, 
we have used (2x1) supercell. Figure [1J a) presents the 
atomic structure of a sample bare armchair Si nanoribbon 
after structural relaxation. Here one can see a (2x1) re- 
construction at the edges, which would be missed if (lxl) 
unit cell was used in the calculation. The edge recon- 
struction is reminiscent of the reconstruction of Si(100)- 
(2x1) surface. In the latter case, two adjacent surface Si 
atoms each having two sp 3 -dangling bonds come closer 
and form a new dimer bond using one sp 3 - dangling bond 
from each atom. At the end, the number of sp 3 dan- 
gling bonds is halved and hence the energy is lowered 
through reconstruction. Similarly, in Fig[TJa) A and B 
atoms come closer to form a bond. Since the nature of 
bonding is modified around B atom, the ABC triangle 
is bowed. At the end, the number of the sp 2 dangling 
bonds is halved. Note that, this kind of reconstruction is 



FIG. 2. (a) Phonon dispersions calculated for hydrogen satu- 
rated ASiNR-9 and for 2D honeycomb structure of Si. States 
that appear above 600 cm -1 are related to Si-H bonds of hy- 
drogen saturated ASiNR-9 and were not shown, (b) Phonon 
densities of states (DOS) of the hydrogen saturated ASiNR-9 
projected to Si atoms at the center and at the edges and also 
to H atoms at the edges. DOS of 2D honeycomb structure of 
Si is also presented for comparison. 



not seen in bare armchair graphene nanoribbons. More- 
over, in graphene and its nanoribbons all atoms lie in the 
same plane, while structures considered here are slightly 
buckled. The separation between adjacent atoms in the 
perpendicular direction to the plane is around 0.4 A and 
0.6 A for Si and Ge nanoribbons, respectively. 

Figure [TJb) presents the atomic structure and bond 
length distribution of a sample hydrogen saturated arm- 
chair silicon nanoribbons. In contrast to bare nanorib- 
bons, saturation by hydrogen lifts the (2x1) reconstruc- 
tion at the edges. In Fig. [TJb) there are n = 9 Si atoms 
forming zigzag chain perpendicular to the nanoribbon 
axis and hence this armchair nanoribbon is classified as 
ASiNR-9. Accordingly the number of Si (or Ge) atoms in 
the primitive unit cell is 2n. Note that the bond length 
distribution is nearly uniform except a sudden decrease 
at the edges. This pattern was also observed in armchair 
graphene nanoribbons. 20 

Figure [2|a) presents the phonon dispersion profile for 
hydrogen saturated ASiNR-9. Also phonon dispersion 
profile of 2D silicon, reproduced from Ref. [14 , was 
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shown for comparison. It is not possible to generate the 
dispersion profile of ASiNR by folding that of 2D Si, but 
the pattern is similar to that found in nanotube disper- 
sion profiles obtained by zonefolding. 28 In the phonon dis- 
persion profile of hydrogen saturated ASiNR-9 all modes 
are real, except for some small imaginary frequencies cal- 
culated for the twisting acoustic mode, TW, near the T 
point. This issue was faced earlier and was attributed to 
the limitations of the computational precision. 29 Thus, 
the structure is predicted to be stable. Computational 
cost of this calculation is very high, so we were not able 
to calculate the phonon dispersions for other ribbons. 
Nevertheless, all ASiNRs have very similar atomic con- 
figuration, and thus they are also expected to be sta- 
ble. The phonon dispersion profile of 2D Ge is similar 
to that of SiP^ But in Ge structure the acoustic and op- 
tic modes are well separated. Also due to softer bonds 
the wavenumbers of Ge structure is halved compared to 
Si. Thus AGeNRs are also expected to be stable, whilst 
exhibiting the mentioned differences. 

Phonon densities of states (DOS) of hydrogen satu- 
rated ASiNR-9 projected to atoms at different locations 
in the nanoribbon are presented in Fig. [5|b). Shown is 
also DOS of the 2D Si honeycomb structure in the same 
figure. DOS projected on Si atoms at the center of the 
nanoribbon is very similar to that of the 2D Si. As the 
width of the nanoribbon increases, this similarity is ex- 
pected to be enhanced. However, DOS projected on Si 
atoms at the edges deviate from that corresponding to 
2D Si. Especially, four optical peaks above 600 cm -1 are 
clearly originating from Si-H bonds at the edges. Also 
modes originating from short Si-Si bonds at the edges 
cause changes in DOS below 600 cm -1 . 



IV. ELECTRONIC STRUCTURE 

In this section, we investigate the electronic structure 
of hydrogen saturated Si and Ge nanoribbons. Nanorib- 
bons having widths from n = 3 to 30 were investigated. 
Figure [3] presents the variation of energy gaps and effec- 
tive masses with the ribbon width given in terms of n. 
Band gaps are direct and located at the T. One can see 
three branches with decaying profiles originating from the 
quantum size effect. Here in ascending order of the band 
gaps, we would like to name these branches with widths 
n = 3/c + 2, 3fc, and 3/c + l as family I, II, and III, respec- 
tively, where k is an integer. This "family behavior" was 
also observed in armchair graphene nanoribbons. 20 This 
trend is explained by foldings of infinite graphene band 
profile, which is easily understood by using a tight bind- 
ing model. The simplest possible model is to assume that 
only the first nearest neighbors interact and have equal 
hopping parameters with self-energies set to zero. But 
this model results in zero band gap for the members of 
the family I. Fortunately, this problem is fixed if a differ- 
ent parameter is used at the edges. 20 As we mentioned in 
the previous section, the Si-Si (and Ge-Ge) bond length 
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FIG. 3. Calculated (a) energy gap and (b) effective mass 
versus ribbon width, n, for hydrogen saturated Si and Ge 
armchair nanoribbons. Filled circles indicate the ab initio 
results while empty circles stand for the results of the tight 
binding fitting. The fitting is performed using only the energy 
gap data. Parameters found from this fitting was used to 
generate the tight binding effective mass data. In each panel 
three branches are observed and named in increasing order of 
band gaps and effective masses as family I, II, and III. 



is apparently smaller at the edges, which implies that 
edge bonds are stronger. That is why, defining a differ- 
ent hopping parameter at the edges reflects the nature 
of the system better. Accordingly, we set the nearest 
neighbor hoping integrals to be t(l + S) at the edges and 
t otherwise. In our model all self energies were set to zero. 
The tight binding parameters were obtained by fitting to 
the first principles results. Results presented in Figj3] (a) 
show that the model used here is successful in reproduc- 
ing the DFT band gap trends of both Si and Ge armchair 
nanoribbons. Since we did not make GW correction to 
the band gaps, the nearest neighbor hopping integrals, 
which are found to be t = 1.03 eV for Si and t = 1.05 eV 
for Ge, are expected to have larger values and should be 
taken as a qualitative result. However, the relative in- 
crease of the hopping integrals at the edges defined by 
S can be taken as a quantitative result. We have found 
S = 0.12 and S — 0.08 for Si and Ge, respectively. Inter- 
estingly, the value found for Si is equal to that of reported 
value for the armchair graphene nanoribbons. To sum up, 
the tight binding parameters determined for Si and Ge 
armchair nanoribbons are t$i = 1-03 eV, tc e — 1-05 eV, 
tsi,ed 9 e = 1.15 eV and t Ge , edge = 1.13 eV. 

Figure [3] (b) presents the effective masses of the first 
conduction band calculated by using the formula 



(1) 



where E(k) is calculated by using both DFT and tight 
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FIG. 4. Band structure and charge density isosurface profiles for several superlattice structures. Solid (green) and dashed (red) 
lines stand for the ab initio and tight binding calculation results, respectively. Structures are labelled, as defined in the text, 
and are given on top of each structure. Top and bottom panels present the results for superlattice structures with AiV = 2 
and 4, respectively. Constituent nanoribbons of super lattices that appear in the same column are members of the same family. 
Charge density isosurfaces calculated by ab initio technique around the T point for two conduction and two valence band edge 
states are ordered in the similar manner as they appear on the energy band structure. 



binding model mentioned above. Note that, parameters 
of the tight binding model are generated by using only 
the band gap information at the T point, but it can re- 
produce the second order momentum derivative and thus 
the effective mass, which is in agreement with that of a b 
initio calculations. The deviation from this agreement 
is seen in the ultimately thin nanoribbons. This is be- 
cause, in these nanoribbons the edges affect the rest of 



the structure and the tight binding model can not be ap- 
plied with the same success. Note that, the effective mass 
trends are similar to the band gap trends and show the 
family behavior. Similar trends in the effective mass were 
also observed in the armchair graphene nanoribbons 
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FIG. 5. Variation of confinements and energy gaps of super- 
lattices with the lengths of constituent parts. Results are 
derived by the tight binding model described in the text. 
Numbers labelling the curves correspond to n\ and 712. The 
lengths, /, of the constituent parts are equal and shown in the 
logarithmic scale. 



V. SUPERLATTICES 

In this section we investigate the electronic structure 
of superlattices formed by periodic junction of ASiNRs 
which have different widths. The variation of energy 
gaps with the ribbon width causes these structures act 
as a multiple quantum well and is expected to lead to 
interesting device applications. In the past, superlattices 
of 2D GaAs-AlAs or Si-Ge heterostructures have been 
extensively studied to realize new generation electronic 
devices. Recently, superlattices of graphene have been 
synthesized experimentally.^ The transport characteris- 
tics and electron confinement have been revealed. 32 Here 
superlattices are formed by the junction of two differ- 
ent nanoribbons. In this study the lengths of constituent 
nanoribbons are taken to be equal. The unequal cases 



were studied for armchair nanoribbons of graphene.Hl We 
let / unit cell of a ribbon with the width of ri\ to make 
a perfect junction with / unit cell of a ribbon with the 
width of 77,2; once joined these structures form a supercell 
having the length of 21. Here ni and ri2 are chosen to be 
odd, so that the ribbons are joined symmetrically. Then 
the resulting structure is relaxed and the lattice constant 
of the whole structure is determined. This superlattice 
structure is labelled by its dimensions as SL(ni,7i2;Z). 
The indices 1 and 2 are arranged in a way that a rib- 
bon with a width rt\ has smaller band gap than a ribbon 
with a width 77*2 . The width difference AN = n\ — ri2 is 
defined for classification purpose. 

Figure [4] presents the band structure and projected 
charge density isosurface plots for sample superlattices 
having AN = 2 and 4 calculated by DFT. The band 
profile is nearly symmetric around the Fermi level and 
one can easily track the bonding and antibonding states 
from the charge density profile. As seen in the charge 
density plots, the band edge states of SL(11,9;4) and 
SL(11,13;4) superlattices are confined in the wide and 
narrow part, respectively. This can be explained by tak- 
ing into account the band gaps of constituent ribbons 
of superlattice structures. Due to the symmetry between 
valence and conduction band edges of the constituent rib- 
bons, the superlattice band line-up is always normal. As 
a result, the part having lower gap acts as a quantum 
well for electrons and holes. That is why, the electrons 
of SL(11, 9; 4) and SL(11, 13; 4) structures are confined in 
the n = 11 part, which for both superlattices is the part 
having the smaller energy gap. 

The situation is also similar for the structures having 
A TV = 4. Electrons of SL(17,13;4) and SL(11,15;4) 
are confined, respectively, in the wide and narrow part, 
which has the lower energy gap. Moreover, the superlat- 
tices constructed by the nanoribbons, which are members 
of the same family, have similar electronic structure. In 
Fig. [4] the superlattices having such common property 
are shown in the same column. Here one can see the 
similarity in the energy band profile of these structures. 
Another way to construct such structures is to increase 
the width of both constituent parts by n = 6k. To ver- 
ify this, we have calculated the electronic structure of 
SL(17, 15; 4) and SL(17, 19; 4), which yielded in the simi- 
lar results as that of SL(11, 9; 4) and SL(11, 13; 4), except 
that in former structures the confinements are less pro- 
nounced because the band gap difference is lower. 

We have to mention that, even though the structures 
SL(9, 13; 4) and SL(15,13;4) have a considerable band 
gap difference between constituent parts, they do not 
have confinement neither in the wide nor in the narrow 
part. This is related to the interface effects, which will 
be discussed in a frame of another simple model. 

Calculations for larger structures with ab initio tech- 
niques are computationally too expensive. So we have 
used the tight binding model mentioned before. In Fig (4] 
one can see that, the tight binding model can reproduce 
the band edge profiles over the whole Brillouin zone. 
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Moreover, the magnitudes of the eigenstates calculated 
by tight binding model (not presented here), mimics the 
projected charge density profile of a given state. To get 
a qualitative picture of how the large structures behave, 
the tight binding model was used. 

Figure [5] presents the energy gap and the band edge 
confinement strength trends generated by the tight bind- 
ing model for the superlattices having n = 9 — 31, 
AN = 2 — 6 and I = 1 — 64. The magnitude squares, 
|^| 2 , of conduction and valence band edge eigenfunc- 
tions generated by the tight binding model are equal. 
The confinement percentage in the narrow (wide) part 
is defined as the sum of the magnitude squares of band 
edge eigenfunctions in the narrow (wide) part multiplied 
by 100 %, noting that the overall sum of the magnitude 
squares is normalized to 1. In Fig. [5] we see that, in 
agreement with the ab initio results presented in Fig. [4| 
the band edge states are confined in the narrow part for 
SL(11, 15; /) and SL(11, 13; /) structures and in the wide 
part for SL(11,9;Z) and SL(17, 13;/) structures. These 
confinements are enhanced as the lengths / of the seg- 
ments increase, because the strength of the well and the 
barrier increase. The trends shown in Fig. [5] also con- 
firm the statement that the superlattices composed of 
the nanoribbons from the same family behave similarly. 
For comparison, in each panel of the Fig. [5] structures 
having such similarity are shown in the same style. Also 
the corresponding family names are presented under the 
figure labels. The similarity can be seen in both confine- 
ment and the energy gap trends. For structures having 
AN — 6, the band edge states are always confined at the 
wide part. This is because, the band gap for structures 
having AN = 6 is lower in the wide part, so the wide 
part acts as a quantum well. 

In Fig. [5j one can see that there is no strong con- 
finement in the wide or narrow part of SL(9, 13; /) and 
SL(15, 13; /) structures. This is also seen in the ab ini- 
tio results presented in Fig. [4j Moreover, the case is 
also true for the superlattices composed of the nanorib- 
bons from the same family. Nevertheless, this can not be 
explained by the band gap difference of the constituent 
nanoribbons, because there are superlattices composed of 
the nanoribbons having similar band gap difference but 
show strong confinement patterns. Actually, plotting the 
linear charge density of former structures along the peri- 
odic direction results in decaying profiles in both narrow 
and wide parts. Thus for these structures, the band edge 
states are localized at the interfaces. This means that, 
the interface acts as a quantum well in these structures. 

To consolidate the effects of the interface, we have to 
make an estimation about its effective energy gap and 
mass. For this purpose, we have chosen the energy gap 
and the effective mass of SL(ni,n2; 1) structure to rep- 
resent the effect of the interface between the nanorib- 
bons with widths n\ and ri2. This information was used 
in construction of a very simple ID model, which qual- 
itatively explains the confinement trends mentioned so 
far. Here the band edge state of a superlattice structure 
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FIG. 6. Model representation of three different superlattice 
structures, (a) Conduction band edge profile when interface 
effects are not included. Here the height of the potential in 
each side is determined by E n , which is the difference be- 
tween the conduction band edge minimum and Fermi level of 
a nanoribbon with a width n. Fermi levels of each structure 
are shown by the dashed lines. The family name correspond- 
ing to each side is written under these lines, (b) Magnitude 
squares (|^| 2 ) of solutions. Dark (blue) and light (green) 
lines represent solutions for the cases where interface effects 
are included and are not included, respectively. The thin line 
joining the panels (a) and (b) corresponds to the geometric 
interface of the superlattices. (c) Conduction band edge pro- 
file when interface effects are included. Here the height of the 
potential at the interface is determined by E ni _ n2 , which is 
the difference between the conduction band edge minimum 
and Fermi level of SL(m,n2; 1) structure. The interfaces act 
as if they are composed of the families I, II, or III as written 
in that region. 



SL(ni, ri2 \ I) is modelled as a single electronic state which 
is under the influence of a periodic potential. The height 
of this potential is taken to be equal to the difference 
between the conduction band minimum and Fermi level 
of each region. The resulting profiles are shown for three 
different superlattice structures in Fig. [6] (a). Here the 
interface effects are not included and the length parame- 
ter of each part is set to 64 times the length of a unit cell. 
Also the mass of electrons in each segment is modelled 
by the effective mass of a nanoribbon with a width n of 
that segment. The magnitude square of the numerical 
solution of this system is presented with a light (green) 
line in the Fig. [6] (b). To include the interface effects, 
the potential height of an arbitrary small region at the 
interface is changed to the difference between conduction 
band minimum and Fermi level of SL(ni, 712; 1) structure. 
The dark (blue) lines in the Fig. [6] (c) represent the new 
effective potential profiles. Solutions of these systems are 
given by dark (blue) lines in the Fig. |6](b). 

Examining the changes induced by inclusion of the in- 
terface effects given in Fig. [6] (b) one can deduce that 
in SL(11,9;64) and SL(11,13;64) structures the inter- 
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face having large effective mass acts as a barrier which 
slightly enhances the confinement strength. In paral- 
lel to this, one can see in Fig. [4] that the charge den- 
sities are more localized in SL(11,9;4) structure than in 
SL(11, 13; 4) structure. Also the enhanced confinement of 
SL(11, 9; I) structure is seen in trends presented in Fig.|5j 
In SL(15, 13; 64) structure, however, the effect of the in- 
terface is dramatic. Since the effective energy gap of the 
interface region is smaller than that of the narrow and 
wide parts, it acts as a quantum well. As a result, we 
have a charge density confined in the interface region 
and decaying in the narrow and wide parts. The decay 
is sharper in the region where the effective potential is 
deeper and the mass is higher. This profile is similar to 
the one found by tight binding model. 



Trends found in this section can be summarized by 
very simple arguments; (i) if the superlattices are formed 
by two different families then the interface acts as if it is 
formed by a family different from both (see Fig. [6] (c)); 
(ii) in this case each part acts differently and electrons 
are confined in the part acting as the family I, which has 
the lowest band gap; (iii) if both sides are composed of 
the same family, then the electrons are confined at the 
wider part; (iv) confinements are enhanced by increasing 
the lengths of the constituent parts. 



VI. CONCLUSION 

We have investigated atomic structure, stability and 
electronic properties of armchair silicon and germanium 
nanoribbons by performing first principles calculations. 
The edges of bare armchair nanoribbons are recon- 
structed. The reconstruction is, however, removed by 
the hydrogen passivation of the edge atoms. It was shown 
that, these nanoribbons exhibit, so called, family behav- 
ior which was explored also in graphene nanoribbons. 
A simple tight binding model proposed for graphene 
nanoribbons was shown to work very well for silicon and 
germanium nanoribbons. Superlattices formed by peri- 
odic modulation of silicon nanoribbon widths were also 
investigated. Modulation of widths in the real space at- 
tributes these structures multiple quantum well proper- 
ties. Specific electronic states are confined in these super- 
lattice structures. Confinements increase with increasing 
the lengths of constituent parts of the superlattice. In 
general, the band edge states are confined in the part of 
superlattice which have the smallest band gap by itself. 
This part can be narrow, wide or the interface part of the 
superlattice. Superlattices composed of the nanoribbons 
which are members of the same family have similar elec- 
tronic structure. Confinement patterns can be explained 
by inclusion of the interface effects in a frame of a simple 
ID effective potential model. It is demonstrated that Si 
and Ge nanoribbons and superlattices constructed there- 
from display features, which may be exploited in future 
nanodevices. 
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